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The Sunyaev— Zel'dovich temperature of the intracluster 
medium 
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ABSTRACT 

The relativistic Sunyaev- Zel'dovich (SZ) effect offers a method, independent of X- 
ray, for measuring the temperature of the intracluster medium (ICM) in the hottest 
systems. Here, using TV-body /hydrodynamic simulations of three galaxy clusters, we 
compare the two quantities for a non-radiative ICM, and for one that is subject both 
to radiative cooling and strong energy feedback from galaxies. Our study has yielded 
two interesting results. Firstly, in all cases, the SZ temperature is hotter than the 
X-ray temperature and is within ten per cent of the virial temperature of the cluster. 
Secondly, the mean SZ temperature is less affected by cooling and feedback than the 
X-ray temperature. Both these results can be explained by the SZ temperature being 
less sensitive to the distribution of cool gas associated with cluster substructure. A 
comparison of the SZ and X-ray temperatures (measured for a sample of hot clusters) 
would therefore yield interesting constraints on the thermodynamic structure of the 
intracluster gas. 
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1 INTRODUCTION 

Currently, nearly all of our knowledge of the intraclus- 
ter medium (ICM) comes from X-ray observations, with 
the latest generation of X-ray satellites, XMM-Newton and 
Chandra, being capable of spatially resolving the density 
and temperature structure independently in low-redshift 
clusters out to tens of per cent of t heir virial radii (e.g . 
Arnaud. Pointecouteau fe Prattll2005l : IVikhlinin et~afl l2005: 



Pratt et a 3 120071 ). Separate density and temperature infor- 



mation is vital for probing the effects of physical processes 
(such as galactic outflows) on the structure of the ICM (e.g. 
from the entropy and pressure profiles), as well as for con- 
structing mass-observable relationships for use in cosmo- 
logical analyses. The latter is currently being driven by X- 
ray s u rveys such as the XM M Cluster Survey ([Romer et al.l 
l200ll : IStanford et HI l200d ) , promising to find many new 
clusters out to redshifts z > 1. 

The ICM is also detectable at centimetre and millime- 
tre wavelengths through the Sunyaev-Zel'dovich (SZ) effect 
(|Sunvaev &; Zel'dovich 1 19721 ) , the Inverse Compton scatter- 
ing of cosmic microwave background (CMB) photons off free 
electrons in the ICM. As the SZ effect is a scattering process, 
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the observed SZ signal is independent of cluster redshift. 
Thus, thousands of new high-redshift clusters are promised 
from upcoming SZ survey s such as that to be performed by 
the South Pole Telescope (jRuhl et al.ll2004h . 

For most clusters, the SZ effect will only measure the 
integrated pressure of the electrons along the line of sight, 
and so its use as a probe of ICM structure will be limited 
without additional X-ray data. The combination of SZ and 
X-ray surface brightness profiles yields a promising new way 
of obtaining cluster temperature information without hav- 
ing to resort to expensive X-ray spectroscopy; such a method 
yields a temperature th at is closer to the mas s- weighted tem- 
perature of the cluster (|Ameglio et al.ll2007h . 

It is possible, however, to measure the temperature 
of the ICM directly in the hottest clusters with multi- 
frequency SZ observations alone, due to significan t rel- 



ativistic effec t s in these sys t ems (e.g. iRephaelil 



Stebbinslll997l: lltoh etai1ll998 



Challi nor &; Lasenbvl 
1998h 



1995 



1998 



Pointecouteau. Giard Sz Barretl I1998T ). In contrast to the 
above, this estimate of the temperature (hereafter refered 
to as the SZ temperature) is weighted by the pressure of 
the electrons along the line of sight. How the SZ and X-ray 
temperature measurements differ systematically, from a the- 
oretical perspective, is the subject of this short paper. Here, 
using iV-body/hydrodynamic simulations of galaxy clusters, 
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we quantify the difference and show that it can be as large 
as a factor of two, primarily due to the dumpiness in the 
ICM, associated with incomplete thermal support. We also 
show that the temperature directly derived from SZ data 
is the more faithful tracer of the underlying cluster poten- 
tial, being within 10 per cent of the virial temperature in all 
clusters studied. 



2 MEASURING CLUSTER TEMPERATURES 

The temperature of the ICM (or any portion of it) is 
most easily measured by fitting a single-temperature plasma 
model to the observed X-ray spectrum. As the emission is 
primarily due to thermal bremsstrahlung (free-free) in hot 
clusters (kT > 2 keV or so), the temperature is determined 
by the energy scale, E ~ kT, above which the emission de- 
cays exponentially. Since a single-temperature model is used 
to fit a multi-temperature distribution, the result is known 
as a spectroscopically-weighted average of the plasma tem- 
perature. 

Recently, iMazzotta et al.1 (|2004 ) calibrated this weight- 
ing using A/"-body/hydrodynamic simulations of galaxy 
clusters and proposed a simple estimator, known as the 
spectroscopic-like temperature 



(i) 



where p is the mass density of a fluid element in volume, 
dV, with temperature, T. Since the weighting factor is w = 
p 2 T~ 3 ^, cooler dense (i.e. lower entropy) gas is weighted 
the most. Naturally, cooling gas increases its density to try 
and re-establish hydrostatic equilibrium, so regions where 
cooling is important (particularly the core) have a strong 
influence on the X-ray temperature of the ICM. 

An alternative, but much harder, measurement of the 
ICM temperature comes from the spectral distortion in- 
duced by a cluster on the CMB, the SZ effect. The fractional 
change in CMB temperature is 
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2.725K (|Mather et al.1 Il999h . T e is the 
electron temperature of the ICM, n e its density, v\ the 
line-of-sight velocity, and dr — n e crTdl the change in op- 
tical depth to Compton scattering along the differential 
line element, dl. The three frequency-dependent pre- factors, 
g(x), /(x, T e ), h(x), where x = hv/kTcMB, allow the separa- 
tion of the SZ signal to be performed (with /(x,0) = for 
all x). 

The first term in equation J2} dominates for reasonable 
assumptions of cluster properties, and is due to the thermal 
energy of the (mainly) non-relativistic electron population. 
The second term is the temperature-dependent correction 
due to relativistic electrons which can be significant for hot 
clusters; it is this term that we will make use of to estimate 
the ICM temperature. Finally, the third term (known as the 
kinetic SZ effect) is due to the bulk motion of the plasma. 

The frequency-independent part of the first term is usu- 
ally expressed as the Compton y parameter 



Table 1. Weights used to define the hot (T > 10 5 K) gas mass- 
weighted temperature, T mw , the X-ray spectroscopic-like temper- 
ature, T s \, and the SZ temperature, Tgz- 
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Following Hansen (2004), we can then define the Compton- 
average of any quantity F as 
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so the Compton- averaged electron temperature is 
1 



(Te) = 
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T e e - dr. 
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The relativistic correction, /(x,T e ), can be approximated 
as f(x,T e ) = T e S(x); higher-order corrections due to ultra- 
relativist ic elect rons are subdominant in a l l but the ver y 
hottest clusters ([Diego. Hansen fe Silkll2003l : lHansenll2004l ). 
In the Rayleigh- J eans limit [x <C 1), g(x) — —2 and 
5(x) = 17k/5m e c 2 (|Challinor &: Lasenbvi[l998h . Also ignor- 
ing the kinetic SZ effect, we can thus approximate equa- 
tion ([2]) as 

AT 



TcMB 



= y[g(x)+6(x) {T e )\. 



(6) 



One can potentially extract (T e ) from mult i- frequency SZ 
data, especially for hot kT > 5keV clusters. This is a differ- 
ent average from the X-ray temperature, weighted instead by 
the gas pressure (so at a given density, the hottest, rather 
than coolest, gas is weighted highest). For simplicity, we 
will denote this temperature Tsz in subsequent discussion. 
A summary of the weights used to define the gas tempera- 
tures is given in Table [1] 



3 CLUSTER SIMULATIONS 

We analyse N- body /hydro dynamic simulations of three 
galaxy clusters, selected from a larger sample already stud- 
ied by Kay et al. (2004; hereafter K2004), to which we 
refer the reader for further details. The clusters were ex- 
tracted from a large cosmological Af-body simulation run by 
the Virgo Consort ium0 and re- simulated a t higher resolu - 
tion and with gas using the GADGET2 code (|Springelll2006h . 
The clusters have similar masses [M v i r = (1.0,0.6,0.9) x 
1O 15 /i _1 M ] and radii [R vir = (2.0, 1.7, 2.0) /i _1 Mpc] at 
z = 0, but were chosen because they have significantly dif- 
ferent merger histories (Powell et al., in preparation). While 
this sample is limited in size (primarily due to the amount 
of CPU time required), it nevertheless gives some indica- 
tion of object-to-object variations on the scale of a typical 
rich cluster. Each cluster was re-simulated twice (see below), 
with 1 — 2 x 10 6 dark matter particles within i? v ir (the dark 
matter and gas particle masses were mjark = 4 x 10 8 /i _1 M© 
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and m g as = 8 x 10 7 H~ 1 Mq respectively). The force reso- 
lution (equivalent Plummer softening length) was fixed at 
e = 10/i _1 kpc in comoving co-ordinates until z = 1, after 
which it was fixed at e = 5/i -1 kpc in physical co-ordinates. 

3.1 Cluster models 

For each cluster, we considered two models for the gas 
physics, taking the total number of simulations performed 
to six. These models are very similar to those studied by 
K2004, so further details may be found there. 

For the first model, labelled non-radiative, the gas was 
subjected to adiabatic forces and an artificial viscosity (to 
generate entropy in shocks) only. This was done using th e 
standard SPH implementation in GADGET2 (|Springelll2006h . 
This model contains the minimum amount of physics re- 
quired to model the formation of the ICM, being driven 
solely by gravitational processes. While a useful baseline, it 
does not match the observational properties of clusters. 

For the second, more realistic, model, labelled feedback, 
the gas was subjected to the following additional processes. 
Firstly, gas particles with temperature T > 10 4 K were able 
to cool radiatively, assuming a metallicity of Z = 0.3Zq. 
Second, gas that had cooled below 1.2 x 10 4 K and reached 
hydrogen densities, 71h > 10 -3 cm -3 , could either form stars 
(i.e. become collisionless) or be reheated to high tempera- 
ture (kT = 25keV), with equal probability (i.e. /heat = 0.5 
in the jargon of K2004). Tests revealed that the choice of 
density threshold and reheating temperature were neces- 
sary for the clusters to have sufficiently-high core entropy 
to lie on the z = X-ray luminosity-temperature relation, 
while the choice of /heat was primarily to produce a sensible 
cooled fraction (on average, only 14 per cent of the baryons 
with in R500 had cooled by z — 0, similar to that observed 
bv lLin. Mohr fc Stanfordlboofo . 

3.2 Estimating the temperatures from the 
simulations 

The various 2D temperature distributions were com- 
puted from maps us i ng t he procedure outlined in 
lOnuora. Kav fe Thomasl (|2Q03h . Briefly, the relevant weight- 
ing was computed for every hot (T > 10 5 K) gas parti- 
cle within a cylinder of length 6i? v ir, and projected radius 
i?vir, centred on the cluster. These weightings were then 
smoothed, and projected along the length of the cylinder 
onto a 2D 800 x 800 pixel array using the projected ver- 
sion of the GADGET2 SPH kernel. Maps of the spectroscopic- 
like temperature, T s i, were computed using the discrete ver- 
sion of equation JTJ) and the SZ temperature, Tsz, using 
equation JSJ. For comparison, we also computed hot gas 
mass- weighted temperature maps (where the weight in equa- 
tion (pQ) is replaced by w = p; see Table [T]). 



4 RESULTS 

We first present various measurements of the mean tem- 
perature within the (projected) virial radius of each simu- 
lated cluster. Besides the X-ray and SZ temperatures, we 
also consider the projected hot gas mass- weighted tem- 
perature, T mw , and the 3D virial temperature, T v i r (see 



Table 2. Various temperature measurements (expressed as an 
energy, kT in keV units) within the virial radius. Column 1 gives 
the label for each simulation; column 2 the virial temperature; 
column 3 the hot gas mass-weighted temperature; column 4 the 
X-ray spectroscopic-like temperature; column 5 the SZ Compton- 
weighted temperature and column 6 the ratio of the SZ to X- 
ray temperatures. Values in brackets denote the ratio of each 
temperature to the virial temperature. 



Cluster 


kT v [ Y 


kT mw 


kT sl 


kT sz 


Tsz/T s i 


NR1 


5.7 


4.2 (0.7) 


4.0 (0.7) 


6.1 (1.1) 


1.5 


NR2 


4.0 


2.7 (0.7) 


2.2 (0.6) 


3.9 (1.0) 


1.8 


NR3 


6.0 


3.8 (0.6) 


2.7 (0.5) 


5.6 (0.9) 


2.1 


FBI 


5.3 


4.4 (0.8) 


5.4 (1.0) 


5.9 (1.1) 


1.1 


FB2 


3.9 


3.3 (0.9) 


3.5 (0.9) 


4.4 (1.1) 


1.2 


FB3 


6.2 


4.5 (0.7) 


4.8 (0.8) 


6.2 (1.0) 


1.3 



iMuanwong et all [2002). the most suitable temperature for 
cosmological applications of clusters. 

The measured temperature^ are given in Table [2] 
The hot gas mass- weighted temperature, T mw , and X-ray 
spectroscopic-like temperature, T s i, are both lower than the 
virial temperature, with the exception of one of the feedback 
clusters (FBI), where T s \ ~ T v i r . In general, the difference 
between these two temperatures and the virial temperature 
is larger for the non-radiative clusters than the feedback clus- 
ters, the most extreme case being NR3, where T s \ ~ T v i r /2. 
Part of this difference can be attributed to projection ef- 
fects, i.e. including gas extending to 3 virial radii from the 
cluster centre, along the line of sight, when calculating T mw 
and T s \. Cooler gas, associated with infalling substructure, is 
more prevalent on the outskirts of clusters, causing a signif- 
icant downward shift in the measured temperatures. When 
the calculation of projected temperatures is restricted to gas 
within i?vir in all directions, the difference between these 
temperatures and T v i r is approximately halved. 

We also show, in Fig. [1] (top panel), how T s \ compares 
to Tvir when the former is measured within a projected ra- 
dius R < i?vir (as is the case with X-ray observations of 
clusters). Thus, the temperature at each value of R is a re- 
sult of integrating over all gas within that radius. At R500 
(~ i?vir/2), the typical outer radius where T s \ can be reliably 
measured with high-quality X-ray data, T s \ is only slightly 
larger (within 10 per cent) than its value at i? v ir- On the 
other hand, the SZ temperature, Tsz, is higher than both 
T mw and T s \ in all clusters studied and is much closer to T v i r 
(within 10 per cent). All clusters show that Tsz converges to 
within 20 per cent of T v i r at R > 0.6i? v ir (~ -ffeoo)- Finally, 
the ratio, Tsz/T s i exhibits very different behaviour between 
the two models, but at small radius in particular, the SZ 
temperature can be significantly (tens of per cent) higher 
than the X-ray temperature. At the virial radius, one non- 
radiative cluster even has Tsz > 2T s i, a massive difference. 

Figure [2] illustrates how the two temperatures (T s i and 

2 In practical terms, the temperatures are at the lower end of 
what might be measurable through the SZ effect, but the overall 
results of this paper should not be significantly affected by using 
hotter clusters. 
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Figure 1. Mean temperature within a projected radius, R/R v { r , 
for the non-radiative (dashed lines) and feedback (solid lines) clus- 
ters. The top panel shows the (X-ray) spectroscopic-like tempera- 
ture profiles, the middle panel the SZ temperature profiles (both 
panels show temperatures in relative to the virial temperature, 
T v i r ) and the bottom panel the ratio of the two temperatures. 
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Figure 2. Azimuthally- averaged temperature profiles for the 
non-radiative (dashed lines) and feedback (solid lines) clusters. 
The top panel shows the (X-ray) spectroscopic-like temperature 
profiles, the middle panel the SZ temperature profiles and the 
bottom panel the ratio of the two temperatures. 
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Figure 4. X-ray (top panel) and SZ (bottom panel) temperature 
distributions (fraction of pixels with a given T/T v i r value, within 
-Rvir) f° r the non-radiative (thin lines) and feedback (thick lines) 
versions of cluster 1. Vertical dashed/dot-dashed lines illustrate 
the mean pixel value for the non-radiative/ feedback models. 



Tsz) vary when measured locally, azimuthally-averaged at 
each radius. Strikingly, the SZ temperature profiles are very 
similar for both models outside ~ 0.4R v ir, suggesting that 
this estimator is insensitive to cluster physics away from the 
core. When the X-ray estimator is used, the feedback clusters 
are systematically hotter than the non-radiative clusters at 
all radii (except in the inner core, where a sharp drop occurs 
due to the presence of gas that has radiated a significant 
amount of its thermal energy). The bottom panel in the 
figure explicitly illustrates the ratio, Tsz/T s i, as a function 
of radius. In both models, the ratio increases with radius, 
although the effect is milder in the feedback clusters, with 
Tsz almost a factor of two higher than T s \ at R v i r . 

The reason why the SZ temperature is hotter than the 
X-ray temperature (and closer to the virial temperature) can 
be seen in Fig. [3] where we present maps of the two quan- 
tities (and their ratio) within i? v ir for the first cluster (NR1 
and FBI). In the non-radiative cluster, the map contains a 
significant amount of dark features, associated with infalling 
cooler, denser gas. Since T s \ weights this gas higher, they 
make a prominent contribution to the temperature of the 
cluster. Note that the nature of the cool gas changes in the 
feedback cluster. There it becomes more spot-like, associated 
with the cores of subhaloes; the combination of cooling and 
strong feedback has removed most of the cool, low-entropy, 
gas in this model. The difference in the amount of cool gas 
in the X-ray and SZ temperature maps is shown clearly in 
Fig. |U where we plot pixel temperature distribution func- 
tions within i?vir- For the SZ distributions, the differences 
at the lowest temperatures are also much less significant. 



5 SUMMARY 

Our results clearly show that the temperature of the ICM, as 
measured through the SZ effect, can differ significantly from 
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Figure 3. Maps of the non-radiative (top panels) and feedback (bottom panels) versions of cluster 1 within their virial radius (2 h~ 1 Mpc). 
From left to right: X-ray spectroscopic-like temperature, Compton- weighted SZ temperature and the ratio of the two. For the first two 
maps, the colour scale illustrates logarithmic temperature values, log(T/K) = [7,8], while for the third, the scale illustrates ratios from 
to 3. 



the X-ray temperature of the cluster, with the SZ tempera- 
ture always being the larger of the two. The X-ray temper- 
ature is weighted by cool dense gas the most, so differences 
between the amount of cool clumpy gas due to, for exam- 
ple, the level of radiative cooling and feedback, are most 
significant. On the other hand, the SZ temperature weights 
the hottest, densest gas (i.e. gas with the highest pressure) 
the most, so is generally a better estimator of the virial 
temperature of the cluster. Consequently, differences due to 
the amount of cooling and heating are less severe. The lower 
sensitivity of the SZ temperature to cluster physics (as com- 
pared to X-ray) is in line with previous comparisons b etween 
X-ray and SZ intensity |da Silva et all l200lL 120041 : iNagai I 



A- ray and intensity ( L< 
120061 : iBonaldi et al.ll2007h 



In practice, the SZ temperature of a cluster is extremely 
difficult to measure, and its applicability limited to the 
hottest clusters, typically kT > 5 keV. It does seem fea- 
sible, however, to design mult i- wavelength experiments that 
are able to constrain SZ tem peratures to ^ 10 — 20 per cen t 
accuracy in the near future (|Knox. Holder fc Church| [2004). 
Whether an instrument will ever be able to measure the 
SZ temperature in enough clusters, and to sufficient accu- 
racy, to be useful for cosmological purpos es (where the viria l 
temperature is the desired quantity, c.f. lEvrard et alJ feoOS) 
remains to be seen, but comparison with the X-ray temper- 



ature of a cluster ought to provide useful information on the 
structure of the ICM. 
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